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f*^) | Abstract 

In this paper, we present the general governing equations of electrodynamics and contin- 
uum mechanics that need to be considered while mathematically modelling the behaviour 
of electromagnetic acoustic transducers (EMATs). We consider the existence of finite de- 
formations for soft materials and the possibility of electric currents, temperature gradients, 
and internal heat generation due to dissipation. Starting with Maxwell's equations of elec- 
tromagnetism and balance laws of nonlinear elasticity, we present the governing equations 
and boundary conditions in incremental form in order to solve wave propagation problems 
■ of boundary value type. 
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1 Introduction 

Electromagnetic acoustic transduction is a technique used to generate and detect mechani- 
cal waves in magnetoelastic conductors for the purpose of nondestructive testing. The de- 
vices, called electromagnetic acoustic transducers (EMATs) or electromagnetic acoustic re- 
ceivers (EMARs), function by exploiting the coupling between electromagnetic and mechanical 
effects in a deformable continuum. The transducer mechanism (as shown in the schematic in 
Figure 1) consists of a current-carrying coil suspended on the specimen under test which is 
magnetized by a large static magnetic field. An AC current is passed through the coil which 
generates a time-varying magnetic field in the specimen. This, on interaction with the static 
magnetic field, leads to an electromagnetic body force that causes the generation of mechanical 
waves in the specimen. In the case of EMARs, the existing mechanical waves in the specimen, 
on interaction with the static magnetic field, produce time- varying magnetic field in the vicinity. 
This generates an AC current in the current-carrying coil which can be measured to analyze the 
mechanical waves in the specimen. Thus, it is very important to understand the exact coupling 
between mechanical and electromagnetic fields for an accurate mathematical modelling of the 
process. 

Several orientations of the specimen and the coil for EMATs have been studied in the recent 
literature. Much of the work (including [lj) is based, among other literature, on the papers 
by Ludwig et al. [2], Ogi [3], and Thompson [1] on the theory and numerical simulation of 
EMATs. Recently, Shapoorabadi et al. [5] have presented the governing equations of EMATs 
using a formulation involving the magnetic vector potential. The list of references is by no 
means exhaustive and there are many more papers and theses on this subject. 
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However, the existing theory of EMATs is, in our view, incomplete in some respects. Many 
theoretical models have been developed and experimental results been obtained for only a si- 
nusoidal steady state of current and displacement (such as [1]). A large volume of literature 
assumes linear elastic deformations in the specimen and a linear coupling between the elec- 
tromagnetic and mechanical phenomena (such as [2], [5]). This is, however, not the case in 
reality with the coupling being nonlinear, as has been presented in the classical texts by Pao 
[6], Eringen and Maugin [7j, and in the recent paper by Maugin [8]. Moreover, with the devel- 
opment of soft polymeric electro- and magnetoelastic smart materials ([9], [TO]) in the recent 
years, one needs to take into account the possibility of finite deformations. In some materials, 
heat generation caused by electric currents may lead to significant changes in magneto-elastic 
properties, hence these effects also need to be taken into account. 

In this paper, we present the equations that govern wave propagation in a finitely deformed 
solid in the presence of electromagnetic fields. On the existing finite deformation, magneti- 
zation, and electric polarization, time dependent small increments in the electric, magnetic, 
temperature, and deformation fields are allowed. We provide a general constitutive formulation 
based on a free energy function that depends on the initial deformation, underlying electromag- 
netic fields, and temperature. This is based on a generalization of the constitutive formulations 
of electroelasticity and magnetoelasticity provided by Dorfmann and Ogden [111 [T2] . The same 
formulation of magnetoelasticity has been used by Destrade and Ogden [13j . and author and 
co-workers [14\ [T5] to study some wave propagation problems in magnetoelasticity. 

In Section 2, we recapitulate the basic governing equations of electromagnetism in continua 
and the laws of thermodynamics as given by Pao [6] and Eringen and Maugin [7]. These are 
used to derive constitutive relations from a free energy function. We then allow for small 
time-dependent increments in Section 4 to study wave propagation. Several 'moduli' tensors 
are also introduced in this section to quantify the couplings in elastic, electric, magnetic, and 
temperature fields. It is expected that the general equations provided here will be solved using 
numerical schemes such as FEM for specific boundary problems related to EMATs. 

2 Basic equations 

The undeformed stress-free reference configuration of a continuous elastic body is denoted by 
B r and its boundary by dB r . Let Bt, the current configuration, be the region occupied by the 
body at time t and dBt its boundary. The material points of body are identified by the position 
vector X in B r which becomes the position vector x in Bt- 

2.1 Kinematics 

The time-dependent deformation (or motion) of the body is described by an invertible mapping 
X that maps points from B r to points in Bt such that x = x(X, t). The function x an d its 
inverse are assumed to be sufficiently regular in space and time. The velocity v and acceleration 
a of a material particle at X are defined by 



where the subscript t following a comma denotes the material time derivative. In this paper; 
Grad, Div, and Curl denote the standard differential operators in the reference configuration 
while grad, div, and curl denote the same in the current configuration. The deformation gradient 
tensor is defined as F = Grad%(X,t) and its determinant is denoted by J = detF, with J > 0. 
Associated with F is the right Cauchy-Green tensor c = F T F. For an incompressible material 
the constraint J = detF = 1, needs to be satisfied. 
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Figure 1: A generic EMAT configuration. Bt and V are the regions occupied by the continuum 
material and a wire of the current-carrying coil, respectively. The corresponding boundaries are 
denoted by dBt and dV. 

2.2 Maxwell's equations 

The well-known Maxwell's equations governing electromagnetic fields in a deformable continua 
(the region Bt in Figure 1) are given by the Dipole-Current Circuit Model of Pao [B] as 

divB = 0, curlE + — = 0, div(e E + P) = p e , (2) 

curl f — B - M ) = J + ^ (e E + P) , (3) 

where B is the magnetic induction vector, E is the electric field vector, M and P are, respec- 
tively, the magnetization and electric polarization of the continuum; J is the electric current 
density, p e is the electric charge density, Eq and po are the electrical pemittivity and magnetic 
permeability of the vacuum respectively. We have used the following field relations in the above 
equations 

D = s E + P, B = (H + M) . (4) 

All the electromagnetic quantities above are defined in the current (Eulerian) configuration. 
They can also be expressed in the reference configuration (Lagrangian form) given by 

B; = JF^B, P, = JF _1 P, E, = F T E, M t = F T M, 

Je = JF' 1 (J - p c v) , pe = Jp e , (5) 

using which the Maxwell's equations can be written in Lagrangian form as 

DivB, = 0, Curl (E, + V x B,) = -B,, t , e Div ( Jc^E,) = p E + DivP;, (6) 

Curl (—cBi - e V x (Jc^Et)) - e^Jc-^j 
V Mo / 

= P,, t + CuTl(M { + VxP { ) + JB- (7) 

In the current-carrying coil of rigid material (region V in Figure 1), the governing equations 
for the electromagnetic fields are 

<9B 1 <9E 
divB = 0, curlE + — = 0, curlB = e e r — + J, divE = 0, (8) 

Ot PoPr Ot 

where the constants e r and p r are the relative electric permittivity and the relative magnetic 
permeability of the current-carrying coil, respectively. 

Outside the material, in a vacuum, the governing equations for electromagnetic fields are 

i9B* f)E* 
divB* = 0, curlE* + — = 0, divE* = 0, curlB* - ^oW)-^- = 0, (9) 
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where we denote a physical quantity in vacuum by a superscript *. 

At the boundary dB r , the following conditions need to be satisfied (see, for example, I15| ) 



N x (E/ + V x Bj - F T E*) = 0, (10) 

N • (Bi - JF _1 B*) = 0, (11) 

N-{s Jc- 1 (E l -F T B*)+P l }=a E , (12) 

N x ( J-Vo 'cB, -M,-Vx (e Jc^E, + P,) - p : F T B*) = - <r £ V s , (13) 



where N is the unit normal to the boundary, V = F v is the particle velocity, V s is the value 
of V at the boundary, a E is the surface electric charge density, and K/ is Lagrangian description 
of the surface current density. We note that F is not defined outside the material and the values 
of F and J calculated on the boundary are used in the above equations. 
The boundary conditions are simpler at &P% and given by 

n x (E - E*) = 0, n (e r E - E*) = 0, 

nx( — -B*)=K, n-(B-B*) = 0. (14) 

2.3 Mechanical balance laws 

Balance of linear momentum, in the absence of a mechanical body force, is given by 

div r + f e = pa, (15) 

where r is the Cauchy stress tensor, p is the mass density, and the electromagnetic body force 
as given in [6j is 

d 

f e = peE + J x B + (gradE) T P + (gradB) T M + — (P x B) + div[v <g> (P x B)]. (16) 

Balance of angular momentum gives 

er + L e = 0, L e = P x E + (M + v x P) x B, (17) 

where e is the third order permutation tensor with components and (er)j = EijhTjk and 
L e is the electromagnetic body couple vector. The above balance equations can be written in 
Lagrangian form using the nominal stress tensor T = JF _1 t as 

DivT + Jf E = p r a, e(FT) + JL E = 0, (18) 
where i E and Ij E are Lagrangian counterparts of the corresponding vectors and are given by 
= J~VeF- t E z + J- 2 (FJ Z ) x (FB,) + F- T [Grad(F- T E z )] T (J- x FP0 
+F" T [Grad(J- 1 FB i )] T (F- T M;) + - [,r 2 (PP,) x (FB/)] 

+J- 1 Div [J- X V {(FP ; ) x (FB/)}] , (19) 

L E = J- 1 (FP,) x (F- T E ; ) + J' 1 (F- T M ei ) x (FB/) . (20) 

On any part of the boundary where the traction is prescribed, the boundary condition may 
be given as 

T T N = t A , (21) 
where t\ is the Lagrangian representation of the traction force. 
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2.4 Energy balance laws and Constitutive relations 

First law of thermodynamics gives the balance of energy as 

dU 



P 



dt 



t : grad v — div q + q + w e 



(22) 



where U is the internal energy, q is the heat flux at the surface, q is the volumetric heat 
generation, the symbol : denotes a scalar product between two second order tensors given in 
component form as r: T = Tj,T,-j, and w e is the electromagnetic power (see, for example, Pao 
[B]) given by 

T „ d /P\ dB 
w e = J e -E e + p— I- • E e - M e ■ — , 

where for a dynamic problem we have defined the effective field variables as 
J e = J-p e v, E e = E + vxB, M e = M + v x P. 
Let $ be the absolute temperature, then we can rewrite (|22p as 

pCp m 

Here c p is the specific heat capacity. On defining the pullback versions of the physical quantities 

q« = JF^q, qi = Jq, w E = Jw e , $i = M, 



q + w e + t: grad v — div q. 



(23) 
(24) 

(25) 



the above equation can be written in Lagrangian form as 

PrCv^ t ( J ~^) = T: Grad(FV) + qi + w E - Div q/ . 
Here, the Lagrangian form of electromagnetic power is given as 



(26) 



(27) 



w E = (FJ ei ) • (F- T E ei ) + Pr 



d /FP, 



dt V P 



- ) + Grad 



FP, 



V 



• (F- T E eZ ) 



-JF~ T M P , ■ 



D_ 

dt 



J _1 FBj) + Grad (J _1 FB ; ) V 



(28) 



If S is the entropy density, then the second law of thermodynamics is given by the Clausius- 
Duhem inequality as 



dS , /q 
p— — h div 



q 



> o. 



dt • — v^y ■& - - (29) 

Substituting equation (|22]) into ([29]) and following the standard Coleman-Noll procedure 
(see, for e.g., [IS]), we arrive at the constitutive relations 



M 



5E eZ ' ~- el dB{ 

where $ is the free energy per unit volume related to U (see, for example, |15j for details) as 



(30) 



(31) 



We note that in the case of an incompressible material, we have an additional constraint 
J = 1 and the constitutive law for stress is modified to 



m d$ , 



(32) 
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p being the Lagrange multiplier associated with the constraint of incompressibility. 

We assume the constitutive laws governing heat flow and electric current flow to be given 
by the Fourier's law and the Ohm's law, respectively, as 



-Kgradi?, J = £E, 



(33) 



where k and £ are positive-definite symmetric second order tensors quantifying the thermal 
conductivity and the electrical conductivity, respectively. We require tensorial values for the 
conductivities since in general for anisotropic materials, the conductivity depends on the direc- 
tion of heat flow or electric current. For isotropic materials, these assume the values k = kI 
and £ = £1. The above equations can be written in Lagrangian form as 



= -JF" 1 KF- T Grad(J- 1 ^) 5 Ji = JF -1 £F- T Ej. 



(34) 



3 Incremental equations 

On the initial motion and underlying electromagnetic fields, we consider an incremental mechan- 
ical motion u(x,t), and increments in electromagnetic fields which are denoted by a superposed 
dot. We emphasize here the departure from convention in using a superposed dot to denote an 
increment rather than a time-derivative for the sake of brevity. 

The incremented forms of the Lagrangian Maxwell's equations (|6j) and (J7J) are given as 



DivB* = 0, CurlE, 



el 



B 



£o Div 



.7(divu)c _1 Ej + Jc^E, - JF~ X (L + L T )F~ T E Z 



PE + DivP;, 



(35) 
(36) 



|U x Curl 



J- 1 {(divu)cBj + 2F T LFB ; + cB z } - e Curl V x (Jc _1 Ej) 
+V x | J(div u)c -1 Ej + Jc^Ei - JF _1 (L + L T )F~ T E,} 
J | J(divu)c" 1 E z + Jc -1 Ej - JF _1 (L + L T )F~ T E^| 



Pi t + CurlM ei + J E . 



(37) 



The incremented Lagrangian quantities are 'pushed forward' to the Eulerian configuration 
and denoted by a subscript '0' after T or l E\ The push-forward relations for the incremented 
fields are given (similar to as 



B 
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J FB 



F~ T M P 



± m = F- T E l p E o = J~ 1 pE, (38) 
using which we can write the Eulerian form of the above incremental equations as 

divB ro = 0, curlE ei0 = [T - (divv)I] B, - B, 0)tl (39) 

e div E = peo + div P /0 , (40) 



/i^curl [{(l + divu)I + 2L}B 



eocurl u t xE-|-vxE — eqE 



curlM ei0 + P,„ t + [(divv)I - r] Pjo+Jeo, 



(41) 
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where 

E = E, + (div u)E - (L + L T )E, (42) 
Edo = F~ T E e , = E, + v x B, + (u,* - Lv) x B, (43) 
j E0 = J^FJe = £E - p EO v - p e (u jt - Lv) . (44) 
The incremented momentum and angular momentum balance equations (|18p are given as 

Div T + J(div u)t E + J?e = p r a, (45) 

£(LFT + FT) + J(divu)L B + JL E = 0. (46) 
where the increments in electromagnetic body force and moment are given by 

h = J" 1 {-(divu)p E F- T E z + p E F~ T E, - />eL t F- t E z + peF~ t E/| 

- J- 1 L T F- 1: [Grad(F~ T E,)] T (FP,) + J _1 F- T [Grad(-L T F- T E, + F~ T E,)1 T (FP,) 
+ J _1 F- T [Grad(F~ T E,)] T [-(divu)FP, + LFP, + FP, 
+ 2J~ 2 (divu)(FJ,) x (FB,) + J~ 2 (LFJ, + FJ,) x (FB,) 

+ J~ 2 (FJ,) x (LFB, + FB/) - L T F~ T [Grad (J _1 FB,)] T (F~ T M,) 

+ F" 



d 



Grad (-J -1 (divu)FB, + J^LFB, + J^FBA (F~ t M,) 
+ F~ T [Grad(J^ 1 FB / )] T (-L T F- T M, + F- T M,) 

2J- 2 (divu) (FP,) x (FB,) + J' 2 (LFP, + FP,) x (FB,) 



+J- 2 (FP,) x LFB, + FB 



J -1 Div 



- J -1 (divu) Div [J~ X V <g> {(FP, x (FB,)}] 
J- X V ® {(FP,) x (FB,)} - J-^div^V {(FP,) x (FB,)} 
' ^LFP, + FP,) x (FB,) + (FP,) x (LFB, + FB,) }] , 

t E = -J-^divu) {(FP,) x (F" T E,) + (F~ T M e ,) x (FB,)} 

+ J' 1 (LFP, + FP,) x (F~ T E,) + J- 1 (FP,) x ^— L T F~ T E, + F _T E, 
+ J" 1 (F~ T M e ,) x (LFB, +FB,) + J" 1 (-L T F" T M e , + F" T M e ,) x (FB, 



(47) 



When updated to Eulerian form, the balance equations (|4"5"]l and (|46|) become 

div T + (divu)f e + f EQ = pu u , 
e(Lr + T ) + (divu)L e + t E0 = 0, 



(48) 

(49) 
(50) 



where f^o and L^o ar e the push-forward forms of the incremental body force and moment, 
respectively, and are given by 



f B0 = -(divu)p e E + p E0 E - p e L E + p c E, - L (gradE) 1 P 



+ 



rad -L A E + Eio 



P + (grad E) 1 - (div u)P + LP + P w 



-2(divu)J x B+ (LJ + J, ) x B + J x ( LB + B m ) — L T (grad B) T M 

T 



+ 



grad ( -(divu)B + LB + B /0 ) M + (gradB) T ( — L 1 M + M zo 
d 



+ 



at i 



2(div u)P xB+ LP + P;o x B + P x LB + B 



(divu) div [v0(Px B)] + div [(v - Lv) <g> (P x B) — (divu)v <g> (P x B) 



+v ® | (l 



LP + P/o ) x B } + P x ( LB + B 



Ho 



(div u) (P x E + M e x B) + ( LP + P l0 ) x E + P x ( — L T E + E 



Wo 



+ M e x ( LB + Bm ) + -L^e + Mrfn x B. 



(51) 



(52) 

The heat equation (|27p can be incremented to give 

Divcu + p rCp ^ - J-^divu)^] =qi+w E + T: Grad(FV) 

+T : Grad(LF V + F V) , (53) 
which when updated to the Eulerian configuration becomes 

div q i0 + P c p-Q t [&lo ~ (div u)i?J = q i0 + w E o + T : grad v + r : grad v. (54) 

We have used the push-forward relations q^o = J -1 Fq/,i?jo = J^&uWeq = J^we, and 
Qio = J~ 1( ii to effect the above transformation. The increment in the electromagnetic power is 
given by 



w E = ( LFJ e/ + F j el ) ■ {F' T E el ) + (Fjy • ( -L T F- T E e/ + F~ T • E el 

LFP Z + FP^ 



+ Pr 



8 / LFPz+FP; 

8t \ p r 

d /FP, 



(F E e j) + Pr 



dt \ p. 



J ( (divu)F" T M eZ - L x F~ x M e ; + F- x Me, 



+ Grad 
+ Grad 

Ttt-Tt 



V + Grad 



FP, 



V 



p r 



V 



-L^F E e ; + F -i • E e 







0_ 

at 



-J- 1 (-(divu)FB ; + LFB Z + FB; 



+Grad ( J _1 FB,) V] - JF- T M d • 

+Grad J -1 |-(divu)FB ; + LFB; + FB;} V + Grad (j^FBj) V 
This is given in Eulerian form as 



(55) 



weo = LJ e + J el0 • E e + J e • — L E e + E 



e/0 



d_ I LP + P;p 

dt \ 



p 



+ grad ( — ) (v — Lv 
P 



^ + P Jt {-p 



| (div u' 



M e — L M e + M, 



• i dB 



M R 



L E e + E e ; 
^{-(divu)B + LB + B ro } 



+ (gradB) (v - Lv)] . (56) 



3.1 Incremental constitutive equations and Moduli tensors 

On incrementing the constitutive equations (|30p . we obtain 



T = AF + BE e i + CB l + m l , Pi = -\^FF + gEa + KBi+X^ij, (57) 

and 

M el = -(K¥ + CE el + MB f + ATA) , (58) 

where we have defined the moduli tensors as 

d 2 $ „ <9 2 $ „ d 2 <& d 2 <5> 
B = — — — , C = — — , £> 



<9F3F' <9E eZ <9F' dB/dF' 

d 2 $ d 2 $ d 2 $ d 2 $ 

d 2 $ d 2 $ _ d 2 $ r _ a 2 $ 

Here A. is a fourth-order tensor, B, C, tF, K, are third-order tensors, </, H, C, AA are second- 
order tensors, and 'D,X,Af are first-order tensors (or vectors). To put them in perspective, for 
a problem involving simpler quasimagnetostatic case as presented in [TT] and [J3], only .A, C, 
and AA. are required. Products in (|57p and (|58p are defined in component form as 

(AF) a i = A a ij3jFjp, (BE e j) a j = B ai \pE e ifi, (C~Bl) a i = C a i\pBlfi, 
(^F)i = F^ajFja, (Q~Ej e i) a = GapEel/3, (H^lja = U-apBip , 

(KF)i = K i]aj F jQ , (CE el ) a = £ Q pE eW , (MB,) a = M a pB^ (60) 

and the following relations hold 

/C = C T , F = B T , C = n T . (61) 

On updating (push-forward to Eulerian configuration) the incremented constitutive equa- 
tions fl57D and ([35]) . we obtain 

T = A)L + BoE ero + C B ro +X> i9 ro , (62) 

P/o = - (BjL + S E e/0 + ^ B ro + X ^o) , (63) 

and 

M eI0 = - (C^L + n^Eeio + MqBio + A/oho) , (64) 
with i?zo = J$i and the updated moduli tensors defined in component form as 

T^oij = FikT>kj, Qoij = J l Fi a FjpQ a p, Hoij = Fi a F^7i a p, 

l 0l = F ik l k , M 0ij = JF^-FfiMafi, N Qi = J- l F^N k . (65) 

Incrementing and updating the constitutive equations (f3l|) . we get 



qzo = — (divu)Kgradi? + 2L/cgradi? — Kgrad i?jo — (div u)i? , (66) 



9 



and 



j,o = £ {(divu)I-L-L T }E + E 
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(67) 



On substituting the incremented updated constitutive equations (|62j) . ([63]) . (f64"|) . ([66]) . and 
([67]) into the incremented updated balance equations (|4"0|) . (|4"T]) . (jMJ), (f50|h and jMJ), we obtain 



Eq div 



E w + {(divu)I- (L + L T )}E 



PEO 



-div ( BqL + C? E e / + ^ B ro + X #. 



/o 



/i 1 curl {(1 + divu)I + 2L} B; - e curl (u ( x E + v x E) - e E 



curl ( C^L + HfjEao + M B l0 + Af ^o J + J 



>E0 



[(div u)I - r] BjL + SoE e/0 + « B W + X ^o 



£jL + ^oEejo + "HoB/o + Zo#io 



div (^L + S EeH) + C B i0 + X>oi?zoJ + (div u)f e + f E0 = pu tt , 
e ( Lr + A)L + B E e/0 + C B l0 + X> ^ ) + (div u)L e + t E0 = 0, 



div 



-(div u)k grad $ + 2Lft grad i? — n grad ( i?jq ~~ (div u)i? 



-pep 
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+ Qzo + ^£0 + T : grad u , 



^o-(divu)^ 
+ [AoI> + B E eZ0 + CoBjo + X>o^o) : grad v, 
along with (|39|) . The governing equations of the incremental fields in the region V are 



divB = 0, e r divE = 0, curlE = - 
while in vacuum, we have 

divB* = 0, div E* = 0, curlE* 
3.2 Incremental boundary conditions 



8B 1 . • 8B 

■—, curlB = J + e e r — , 

ot HoHr ot 



dB* dB 
-, curlB = EofiQ 



dt 



ot 



(68) 



(69) 

(70) 
(71) 



(72) 



(73) 



(74) 



At the boundary dB r of the continuum and vacuum, the incremental form of the boundary 
conditions (|10p - (|13|) is given as 



Nx (Ej + VxBi+VxBj - F r LE* - F T E* 



0, 



N 



• (Bi - J(divu)F _1 B* - JF^L^B* - JF _1 B 



B* = 0, 



(75) 
(76) 



N • {eoJc- 1 (±i - F T L T E* - F T E*) + e ^(div u)^ 1 (E, - F T E*) 

+e J (F- 1 L- 1 F- T + F- 1 l- T F- T ) (E, - F T E*) + P/} = a E , (77) 
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N x | J~Vo X cB ; - J _1 (divu)Ato 1 cB l + J~Vo 1 (F Tlf + F T L T F) B, 

—Ml - V x (e Jc -1 Ej + P^) - e V x ^J(divu)c _1 E, + Jc^E, 

+J (F" 1 L" 1 F- T + F _1 L _T F _T ) Ei + P,) - /i " 1 F T L T B* - Mo 1 F T B*} 

= K l -& E V s -a E V s . (78) 

We recollect the Nanson's formula n da = JF~ T N dA connecting reference and current 
area elements dA and da, where n is the unit outward normal to dBt- The surface current 
and electric charge densities in the two configurations are related by K/ = F~ 1 ~Kda/dA and 
oe = (J e da/dA. Using these relations and the relations (|38|) . we rewrite the above incremental 
boundary conditions in the current configuration on dBt as 

n x (E l0 + (v - Lv) x B + v x B l0 - LE* - E*) = 0, (79) 
n • (B l0 - (div u)B* - L X B* - B*) = 0, (80) 

n • {e (e, - L T E* - E*) + e Q (div u) (E - E*) + P l0 

+e {L- 1 + L- T )(B-B*) = a E0 , (81) 

n x {/z^Bio - Mo ^divujB + Mo 1 ( L + L T ) B - M m - (v - Lv) x (e E + P) 
-e v x ((divu)E + E /0 + (L- 1 + L~ T ) E + P, ) - Mo ^ T B* - m *B*} 

= K w + L- 1 K + ( j e Lv s , (82) 

where we have defined K^o = EK/ dA/da and &eo = &E dA/da. 

At the boundary dV, incremental form of the boundary conditions (|14p is given as 

0, n • (e r E - E*) = 0, 

0, n • (B - B*) = 0. (83) 

Thus, a generic EMAT problem requires the solution of equations ([39]) and ([68]) - ([72]) in Bt, 
equations (f73|) in V, and equations ([71]) in vacuum using the boundary conditions ([79]) - ([82]) at 
dB t and (p]) at dV. 

4 Concluding remarks 

This paper reviews the existing theory of electromagnetic interactions with a solid continuum. 
In particular, we analyze the problem of wave propagation in a finitely deformed elastic solid 
with an underlying electric, magnetic, and temperature field. Unlike the existing literature, we 
generalize a nonlinear theory of elasticity to include electromagnetic effects. On a finite defor- 
mation, magnetic, and electric field, the equations are linearized to consider wave propagation 
and several moduli tensors are introduced. In addition to filling a gap in literature, it is ex- 
pected that this work will lead to an accurate mathematical modelling of the 'Electromagnetic 
acoustic transducers' in particular. 



n x (E - E*j = 
nx (*-B*) = 

V Mr I 
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